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1.  Objectives  and  Accomplishments 

The  main  objective  of  this  program  that  was  funded  in  conjunction  with  the  parent 
program  (AFOSR  No.  F49620-96-1-0309)  was  the  feasibility  study  for  modulating  the 
refiactive  index  in  HI-V  semiconductor  multiple-quantum-well  (MQW)  materials. 

This  program  resulted  in  one  Ph.D.  Thesis  authored  by  M.  Hansen,  entitled  "High 
intensity  phase  modulation  using  absorption  saturated  multiple  quantum  wells,"  Ph.D. 
Dissertation,  Electrical  Engineering  (Applied  Physics),  UCSD. 

Traditionally,  thin  film  IH-V  devices  have  been  )hi*iied  to  the  use  of  absorption 
modulation  at  low  optical  intensity.  Although,  the  index  modulation  in  MQW  material  is 
considerably  higher  than  in  other  electro-optical  materials,  this  effect  is  masked  by  the 
absorption  leading  to  phase  modulated  devices  with  very  high  insertion  loss. 

During  the  course  of  this  program,  we  have  developed  a  measurement  setup  to  measure 
index  of  refraction  in  DI-V  thin  film  semiconductors  and  investigated  theoretically  and 
experimentally  the  refractive  index  modulation  in  MQW  and  DX  materials.  We  have 
shown  that  although  small  in  magnitude  index  modulation  exi^s  above  absorption 
saturation.  Although  we  spent  considerable  effort  in  growing  DX  niaLterials  we  were  not 
able  to  demonstrate  conclusively  optical  modulation  in  this  type  of  miaterials 

In  this  report,  we  summarize  the  the  refi:active  index  modulation  in  MQW 
materials  beyond  absorption  saturation,  describe  the  phase-shift  interferometry  method, 
and  detail  the  measurement  results.  We  also  expand  on  our  research  on  DX  center 
materials. 

2.  Theoretical  Background 

The  refractive  index  of  a  material  is  the  macroscopic  representation  of  light  matter 
interaction  when  a  photon  field  is  passing  through  a  medium.  The  interaction  takes  a 
form  of  photon  absorption  and  re-emission  by  atomic  dipole  oscillators.  When  there  is  no 
energy  loss  in  the  passing  photon  field,  the  net  result  of  the  interaction  is  a  phase  delay 
represented  by  the  refi'active  index.  A  permanent  energy  ’oss  occurs  when  oscillators  are 
excited  to  any  higher  energy  states  resulting  fi'om  the  absorption  of  photons.  This  usually 
occurs  at  the  oscillator  resonant  energy.  The  refractive  index  is  larger  near  the  resonant 
energy  as  well  because  the  light  matter  interaction  is  strongest  at  this  energy. 

The  electron  energy  states  are  quantized  in  MQW  material,  and  the  resonant 
conditions  occur  at  each  quantized  levels.  The  exciton  (spatially  bounded  electron/hole 
pair)  resonant  energies  are  slightly  below  each  quantized  energy  level.  Corresponding 
refractive  indices  are  higher  near  the  exciton  energy  levels  and  photon  absorption  occurs 
at  these  levels.  Thus  there  are  specific  energy  levels  corresponding  to  exciton  oscillators, 
for  which  absorption  and  refraction  in  a  MQW  material  can  be  associated  with. 

Absorption  and  index  are  coupled  in  an  absorbing  material  through  the  Kronig- 
Kramers  relation.  Since  the  absorption  coefficient  in  MQW  material  can  be  modulated 
electrically  through  the  Quantum-Confined  Stark  Effect,  the  density  of  excitons  is  thus  a 
function  of  applied  electric  field.  Since  the  density  of  excitons  also  governs  the  refractive 
behavior  of  the  material,  the  refractive  index  is  also  modulated  by  the  electric  field. 
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Wavelength  (A) 

Figurel:  Lx)w  optical  intensity  electroabsorption  spectrum  for  InGaAs/InAlGaAs  MQWs. 

Notice  the  shift  of  the  exciton  absorption  peak  to  lon^gr  wavelength  with 
applied  field.  Broadening  and  reduction  of  absorption  p^  height  also  occur 
with  applied  field  due  to  the  reduction  in  electron-hole  wavefunction  overlap 


The  absorption  coefficient  spectra  can  be  measured  in  function  of  applied-electric-field 
conditions.  Figme  1  shows  measured  results  in  an  InGaAs/InAlGaAs  MQW  material. 
Refiractive  indices  can  be  calculated  following  the  classical  Kramers-Kronig  relation. 
Figure  2  shows  the  calculated  reflective  index  fi'om  the  0-V  curve  in  Figure  1.  Similarly, 
one  can  repeat  the  calculation  for  each  curve  in  figure  1  and  obtain  the  reflective  index 
spectrum  at  each  field  condition.  Figure  3  shows  the  calculated  refi-active-index  spectra. 
The  index  modulation  can  thus  be  obtained  with  the  applied  electric  field  at  a  specific 
optical  wavelength.  For  example,  the  expected  index  modulation  is  -0.01  at  the  .8tV 
applied  voltage  at  the  wavelength  1.064pm. 
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Figure  2;  The  measured  exciton  absorption  resonance  was  isolated  from  the  absorption 
spectrum  (solid  line),  then  its  effect  on  the  index  as  an  isolated  oscillator  was 
estimated  via  a  Kramers-Kronig  transform  (dashed  line). 


Wavelength  (A) 

Figure  3:  Electro-refraction  spectrum  estimated  from -the  electro¬ 
absorption  spectrum  through  the  use  of  a  Kramers-Kronig 
transform 

However,  in  MQW  materials,  the  exciton  absorption  can  saturate  at  rather  low 
intensities.  It  is  then  interesting  to  consider  the  refractive  index  of  a  MQW  material 
where  the  exciton  absorption  is  saturated.  After  the  absorption  is  saturated,  the  material 
behaves  more  like  a  dielectric  material  with  a  constant  refractive  index.  The  value  of  this 
constant  depends  on  the  nature  and  the  density  of  the  oscillators.  The  density  of  exciton 
oscillators  should  be  at  its  maximum  when  absorption -is  saturated,  and  the  maximum 
value  is  a  function  of  the  absorption  coefficient  before  the  saturation.  ~ 

Several  phenomena  occur  when  the  absorption  saturates  in  a  material  as  a  result  of 
high  concentration  of  photo-generated  particles.  They  affect  the  absorption  coefficient, 
the  refractive  index,  and/or  the  internal  electric  field.  These  effects  include  Bandfilling, 
Bandgap  Renormalization,  Thermal  Effects,  and  Electric  field  screening  and  can 
significantly  affect  the  measurement  of  the  refractive  index  imder  saturated  absortion 
condition  in  MQW  materials.  These  effects  have  been  studied  in  detail  in  the  parent 
program. 

3.  Index  of  refraction  measurement 

To  measure  the  refractive  index  modulation  above  saturation  intensity,  a  method  of 
phase  shift  interferometry  was  used.  To  avoid  the  thermal  effect,  a  pulse  laser  with 
picosecond  pulse  duration  was  used  in  the  experiment.  The  interference  fiinges  were 
generated  by  a  wedged  sample  of  surface-normal  reflection-mode  modulator.  Interference 
fnnges  in  the  optical  window  of  the  MQW  modulator  were  magnified,  imaged  onto  a 
CCD,  and  captured  by  a  computer.  Then,  the  phase  of  the  fnnge  pattern  was  determined 
row  by  row  in  the  image  through  Fourier  analysis.  Fringe  phases  were  compared 
between  pairs  of  images  captured  with  the  modulator  in  the  biased  and  unbiased  states. 


The  change  in  the  index  of  refraction  between  these  two  states  could  then  be  calculated 
from  the  measured  phase  change. 

3.1  Sample  Preparation 


The  material  used  in  these  experiments  was  designed  to  operate  at  wavelengths 
compatible  with  commonly  available  high  power  lasers,  such  as  YAG  and  YLF  lasers, 
near  1.0fi4pm.  InGaAs/InAlGaAs  quantum  well  material  can  be  grown  on  GaAs  wafers 
to  show  an  exciton  resonance  at  the  desired  wavelength.  For  operation  at  or  near 
1.06pm,  however,  the  indium  concentration  in  the  hiGaAs  well  material  must  be  near 
20%.  This  represents  close  to  a  2%  lattice  mismatch  with  respect  to  the  GaAs  substrate. 
Strain  relief  occurs  in  the  form  of  dislocations  that  propagate  to  the  surface,  degrading 
device  performance.  Buffer  layers  can  mitigate  strain  effects  caused  by  the  lattice 
mismatch.  Such  layers  are  grown  underneath  the  quantum  well  material  to  gradually 
change  the  lattice  constant  to  that  Of  the  well  material.  As  has  been  demonstrated,  the 
techniques  of  using  linear  graded,  step  graded,  and  superlattice  buffers  have  been 
effective  in  filtering  the  threading  dislocations  generated  from  the  lattice  mismatch.’^’’^ 
The  buffer  layers  prevent  a  majority  of  the  threading  dislocations  ^m  propagating  up 
through  the  material  and  into  the  active  device  layers. 


PIN  diode  modulators  were  fabricated  in  the  epitaxial  layers  on  the  front  surface  of 
the  GaAs  wafer  sample  using  common  lithographic  and  wet  etching  techniques.  The 
backside  of  the  sample  was  mechanically  polished  with  an  intentional  wedge  to  generate 
interference  fringes  at  a  particular  spatial  frequency.  Then  a  high  reflectivity  (at  near-IR 
wavelengths)  gold  mirror  was  thermally  evaporated  onto  the  backside  of  the  sample.  The. 
result  is  a  reflective  mode  surface  normal  MQW  modulator  with  an  internal  cavity  to 
generate  interference  fringes. 


Optical  Window 


n-Contact 


p-Buffer  Layer 


p-Contact 


GaAs  Substrate 


Figure  4:  Cross-section  view  of  the  PIN  modulator 


The  MQW  material  was  processed  into  PIN  ring  diodes,  with  an  optical  window  in 
the  center  of  a  ring  contact.  See  Figure  4  for  a  diagram  of  the  diode  structure.  Wet  etch 
mesa  isolation  was  used  to  electrically  separate  the  modulators.  A  common  ground 
plane,  the  p-doped  buffer  layers,  was  used  for  one  electrical  contact.  The  other  contact. 


Figure  5:  Several  sizes  of  ring  diodes  fabricated  on  the  InGaAs/InAlGaAs  MQW 
material.  Coherent  illumination  readily  shows  the  straight  interference 
fringes  generated  in  diode  windows  due  to  the  sample  wedge  angle.  An 
electrical  probe  tip  is  visible  at  the  lower  right.  The  straight  bar  leading  into 
the  large  size  modulator  is  a  break  in  the  ground  plane  metallization  caused 
by  shadowing  of  the  sample  holder  clip  during  the  metal  evaporation. 
Interference  fringes  can  be  seen  in  this  area  as  well  as  the  diode  windows 


the  n-doped  cap  layer,  was  located  on  top  of  the  mesa.  A  typical  lithographic  and  wet 
etching  procedure  was  used  to  fabricate  the  devices.  Electrical  contacts  were  made  using 
annealed  AuGe  and  CrAuZn  metallization. 

The  ring  diodes  were  fabricated  in  several  sizes  on  the  sample,  as  can  be  seen  in 
Figure  4.  This  image  was  formed  with  coherent  illumination  at  1064nm.  The  straight 
interference  fringes  due  to  the  wedged  sample  are  easily  seen  in  the  ring  diode  windows. 
The  front  surface  of  the  modulator  has  a  reflectivity  of  roughly  30%  from  the  air- 
semiconductor  interface,  which  interferes  with  the  light  returning  from  the  back  reflector. 
The  wedge  between  the  two  surfaces  generates  straight  interference  fringes  when  the 
device  is  illuminated  with  coherent  light. 
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The  fringe  phase  determination  method,  which  is  described  in  detail  below,  requires 
approximately  10  fringes  to  obtain  an  accurate  measurement’^.  Therefore  an  optical 
window- diameter  of  275pm,  was  chosen  as  a  target  size  for  the  10  fringes. 

In  order  to  ensure  a  fringe  spacing  of  27.5  pm  to  fulfill  the  10-fiinge  requirement,  the 
backside  of  the  sample  was  polished  with  a  wedge  angle  of  approximately  0.6®.  The 
wedge  angle  can  be  derived  from  the  desired  fiinge  spacing  (10  Mnges  per  275pm)  and 
the  phase  change  for  one  complete  fringe  (27i).  The  wedge  angle  can  be  found  using  the 
relation 

-ife)  (1) 

where  A,f  =  Xp/n,  Xp  is  the  free  space  wavelength  of  light  and  n  is  the  index  of  refraction  of 
the  wedge  material,  x  is  the  total  length  of  the  fringes,  N  is  the  number  of  fringes,  and  0 
is  the  wedge  angle  as  in  Figure  6. 

The  accuracy  of  the  wedge 
polishing  process  was  optically 
monitored  with  the  sample  mounted 
on  the  polishing  chuck.  Adjustments 
to  the  process  could  then  be  made 
with  a  minimum  of  handling  of  the 
device.  After  the  polishing  procedure 
was  completed,  a  high  reflectivity 
gold  mirror  was  thermally  evaporated  Figure  6:  Geometry  of  an  optical  wedge 
onto  the  back  surface.  Approximately 
2000A  of  gold  was  evaporated  onto 
50A  of  chromium  used  as  a  sticking  layer. 

3.2  Measurements 

Two  sets  of  experiments  were  taken.  First,  photocurrent  measurements  were 
performed  to  estimate  the  saturation  intensity  for  the  MQW  diodes.  Then,  a  high 
precision  phase  shift  measurement  were  made  to  directly  measure  the  change  in  index  of 
refr:uction  in  the  MQW  layers  under  high  intensity  pulsed  operation. 

2.2.1  Photocurrent  Measurements 

The  measurements  were  performed  under  pulsed  optical  input  to  determine  the 
saturation  power  of  the  MQW  modulator  under  different  bias  conditions.  The 
photocurrent  was  monitored  as  a  function  of  input  optical  intensity  and  applied  bias.  The 
low  intensity  linear  unsaturated  regime  and  high  intensity  saturated  regimes  are  clearly 
delineated. 

Using  the  largest  size  diode  on  the  MQW  sample,  window  area  0.484  mm^  (see 
Figure  4),  the  current  generated  from  the  lOHz,  35ps  laser  pulses  was  measured  with  a 
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high-resolution  ammeter.  A  Hewlett-Packard  Model  3457A  digital  multimeter  with 
nanoamp  resolution  was  used.  The  DC  current  setting  of  the  instrument  was  \ised,  and 
the  photocurrent  generated  in  the  reverse  biased  diodes  from  the  short  laser  pulses  was 
enough  to  register  on  the  ammeter  on  its  most  sensitive  (nA)  scale.  By  using  the  DC 
setting,  the  value  measured  is  the  integrated  current  over  the  sampling  period  of  the 
ammeter.  This  approach  does  not  quantify  the  high-speed  characteristics  of  the  current 
pulse,  however  it  does  provide  auseftil  measure  of  photocurrent  saturation. 

In  order  to  estimate  the  saturation  intensity  for  the  measured  photocurrent  curves,  the 
data  was  fitted  with  a  Lorentzian  saturation  equation,  1/(1+I/Is),  in  the  form  of 


abs 


(K)  = 


IikiV) 


1  + 


isiy) 


(2) 


where  7/  is  the  incident  optical  intensity  and  /,  is  the  saturation  intensity.  The  k(V) 
parameter  is  related  to  the  absorption  coefficient  at  low  opticrl  intensity  oto,  the 
absorption  length  of  the  modulator  material  L,  and  the  modulation  voltage  V  through 
k(V)  =  l-exp[-ao(V)L].’^ 

The  measured  data  is  shown  in  Figure  7  along  with  fits  to  the  experimental  data  using 
Equation  (2). 


Figure  7:  Photocurrent  saturation  characteristics  with  pulsed  laser  input 


Table  (1)  Saturation  intensities  for  optical  excitation  with  35ps  laser  pulses 


2.5 

.088 

s'oxlO^ 

5 

.173 

4.4x10^ 

7.5 

.407 

4.7x10=’ 

10 

.501 

7.0x10^ 

12 

.653 

7.9x10^ 

Table  1  Fists  values  for  /,  found  from  the  curve  fits  to  the  experimental  data. 
Saturation-intensity  increasing  with  applied  bias  is  evident.  This  suggests  faster 
carrier  sweepout  from  lowered  barriers  from  the  electric  field.  Values  of  /,  from 
100  to  632  kW/cm^  have  been  observed  by  several  sources  for  similar 
InGaAs/GaAs  MQW  materials. The  results  presaited  here  are  in  line  with  those 
reported  values.  These  values  are  considered  in  the  high  range  for  MQW  absorption 
modulators.  High  saturation  intensities  can  be  attributed  to  fast  carrier  escape  rates 
and/or  fast  recombination  times.  MQW  materials  with  a  high  density  of 
dislocations  from  strain  relaxation  tend  to  have  low  quantum  efficiencies  due  to 
carrier  recombination  at  dislocation  sites  in  the  wells.  This  fast  recombination 
decreases  the  effective  carrier  lifetime,  so  modulators  with  hi^er  dislocation 
densities  tend  to  have  higher  saturation  intensities. 

3.2.2  Interference  Fringe  Shift  Measurement 

To  measure  the  change  in  index  of  refraction  in  the  modulator  we  used  fringe  shift 
interferometry  method.  Fringes  result  from  the  reflections  of  a  laser  beam  from  the  front 
and  back  surfaces  of  the  modulator.  The  shift  of  these  fiinges  that  result  from  an  applied 
electric  field  can  reflect  a  change  in  index  of  refraction.  The  fringes  are  imaged  onto  a 
CCD  and  are  captured  with  and  without  an  applied  electrical  bias  on.  the  modulator.  The 
phase  shift  of  the  interference  fringes  between  the  biased  and  imbiased  conditions  is  then 
measured  with  a  Fourier  analysis  of  each  fiinge  pattern.  The  change  in  index  of 
refiraction  in  the  modulator  can  then  be  determined  once  the  phase  shift  is  known.  Several 
averaging  techniques  as  well  as  a  simple  statistical  analysis  method  were  employed  to 
reduce  noise  and  give  an  estimate  of  the  significance  of  the  results. 

Allo^ving  for  averaging  and  statistical  analysis  techniques,  the  measurement 
technique  provides  for  a  minimum  detectable  phase  shift  on  the  order  of  X/1000  with  a 
sample  size  of  50,  and  X/2500  or  better  with  larger  sample  sizes. 

This  technique  for  direct  measurement  of  the  change  in  the  index  gives  the  sign  as 
well  as  the  magnitude  of  the  change.  It  thus  yields  significant  clues  to  the  physics 
supporting  the  experiment.  This  hopefully  validates  and  simplifies  the  interpretation  of 
the  experimental  results. 

The  optical  setup  used  in  these  experiments  was  designed  with  several 
requirements  in  mind.  The  imaging  system  needed  to  deliver  variable  optical  power  to 
the  diode  modulator.  The  modulators  were  electrically  probed,  so  a  variable 
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magnification  system  was  necessary.  The  sample  is  aligned  into  the  system  and 
electrically  probed  under  low  magnification.  The  magmfication  is  then  increased  so  that 
the  fiinge  pattern  in  the  diode  window  completely  fills  the  CCD,  giving  the  highest  pixel 
count  per  complete  fiinge.  A  high  degree  of  mechamcal  stability  was  also  required  since 
the  fiinge  shift  measurement  is  sensitive  to  vibrations. 

The  laser  used  in  these  experiments  was  a  Continuum  PY61C  pulsed  YAG  laser.  It 
genera^cJ  35ps  pulses  at  a  lOHz  repetition  rate.  The  laser  pulse  was-measured  to  have  a 
top  hat  intensity  profile.  Thus  there  is  a  relatively  constant  intensity  across  the  laser  spot, 
as  compared  to  a  typical  Gaussian  intensity  profile. 


Figure  8:  Layout  of  optical  setup  used  in  the  precision  fiinge  shift  measurement 


The  camera  used  to  capture  the  fiinge  images  was  a  Princeton  Instruments  Model 
TE/CCD-1317-K/1.  This  16  bit  camera  uses  a  cooled  CCD,  operating  temperatoe  of 
10°C,  with  1316x1034  pixel  resolution.  The  camera  is  controlled  via  a  Macintosh 
computer  and  IPLabSpectrum  image  processing  software  from  Signal  Analytics. 

The  MQW  modulator  sample  was  mounted  on  a  six-axis  positioner:  XYZ,  pitch, 
yaw,  and  rotation.  The  MQW  sample  was  probed  electrically  using  micro-positioners 
and  standard  25  pm  probe  tips.  The  probes  were  connected  via  50  ohm  coaxial  cable  to  a 
HP  E3631A  DC  variable  power  supply  and  a  HP  3457A  current  meter.  The  current 
meter  was  used  in  this  experiment  to  verify  electrical  contact  had  been  made  with  the 
probes  and  to  monitor  the  I-V  characteristics  of  the  diode  before  and  after  the  fiinge  shift 

experiment.  ' 

During  the  measurement,  the  modulator  sample  is  flood  illuminated  with  35ps 
lOHz  YAG  laser  pulses.  The  laser  used  could  provide  ample  optical  power  in  order  to 
operate  the  MQW  modulator  well  beyond  its  continuous  wave  (CW)  absorption 
saturation  point.  The  YAG  laser  was  used  at  an  average  power  output  of  300mW.  With 
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a  measured  laser  spot  size  of  13.2mm  in  diameter,  this  translates  into  a  pulse  intensity  of 
6x10*  W/cm^.  A  typical  CW  saturation  intensity  for  a  MQW  modulator  is  on  the  order  of 
10^  W/cm^,  ”  so  there  is  ample  optical  power  available  fix)m  this  laser,  even  when 
unfocused.  The  saturation  limit  was  determined  through  monitoring  of  the  photocurrent 
as  was  described  above  in  the  section  related  to  photocurrent  measurements. 

Each  image  captured  in  these  experiments  r^resents  a  650ms  exposure.  Thus  there 
are  six  optical  pulses  integrated  into  each  image,  given  the  laser’s  lOHz  rep  rate.  This 
multiple  pulse  integration  technique  was  used  to  reduce  the  effects  of  any  intensity 
variation  from  pulse  to  pulse. 

Images  were  captmed  in  groups  of  50  via  automated  image  acquisition.  Each  set  of 
50  images  took  approximately  3  minutes  to  gather.  This  automated  approach  was  taken 
to  speed  data  collection  in  order  to  reduce  error  from  component  drift.  Special  effort  was 
taken  to  minimize  shifting  of  the  image  of  the  fringes  on  the  CCD  due  to  reasons  other 
than  a  change  in  the  index  of  refraction  of  the  sampler.  The  use  of  a  fast  data  collection 
technique  proved  to  be  critical  in  bypassing  slowly  varying  drift  components  of  error.  To 
accurately  determine  the  spatial  frequency  of  the  fringes,  a  one-dimensional  FFT  was 
tperformed  for  each  row  of  a  single  image.  Thus  100  FFTs  were  performed.  As  a  noise 
reduction  technique  the  coliunns  of  the  FFT  data  wwe  summed.  So.  a  one  dimensional 
FFT  spectrum  was  obtained  from  the  summation  of  the  100  FFTs-frxnn  the  entire  image. 

To  minimize  error  in  the  phase  calculationdetermination  of  the  spatial  frequency 
should  be  performed  to  an  accuracy  of  ±0.1  or  better.^’  The  FFT  function  returns  values 
at  integer  frequency  intervals,  so  an  interpolation  is  performed  to  find  the  frequency  of 
the  peak  in  the  FFT  data  when  it  falls  between  integer  values.  In  this  manner,  the  spatial 
frequency  of  the  interference  fringe  pattern  is  easily  determined  to  within  O.I. 


A  profile  of  the  interference  fringes  generated  in  the  window  of  the  MQW  modulator 
can  be  written  as 

Ii(x)  =  W(x)[l+cos(27ifx-4,)]  (3) 

where  x  is  the  spatial  coordinate,  and  f  and  <j>i  are  the  spatial  frequency  and  intrinsic  phase 
of  the  fringes.  The  window  function  of  the  fringe  pattern  truncating  the  signal  is 
represented  by  W(x).  The  fringe  profile  after  applying  a  bias  voltage  (assuming  there  is 
an  index  change  in  the  well  material)  can  be  expressed  as 

l2(x)  =  W(x)[l  +  cos(27ifx  -  <{>1  -  5(|>)]  (4) 

where  the  change  in  phase  is  represented  by  frcj). 

To  find  the  phases  of  the  interference  ^nges,  the  Fourier  cosine  and  sine  integrals  are 
calculated  as 

I 

C,  =  ^I.{x)cos{l7rfx)dx 


(5) 


it 


Si  =  jli(x)smi2f^)dx 
0 

where  the  range  of  the  hinges  is  normalized  and  i=l,2.  From  the  Fourier  sine  and  cosine 
functions  the  phase  of  the  fringe  pattern  can  be  foimd  with  the  arctangent  function: 

<l)i  =  tan-’(Si/Ci)  (6) 

where  <j>2  is  actually  (|>i+5^.  Therefore,  the  phase  change  can  be  foimd  as  the  difference 
between  the  phases  (|»i  and 

There  are  100  rows  of  fringe  intensity  data  in  each  image  captured  by  the  computer. 
The  phase  of  each  row  is  found  using  the  discrete  version  of  Equation  (6): 


where  N  is  the  number  of  pixels  per  row  and  the  quantity  27f/N  is  the  normalized  spatial 
frequency. 

Nakadate*^  used  a  similar  method  of  analysis  to  measure  phase  shifts  in  a  common  . 
path  interferometer  which  used  Young’s  fringe  generation  from  a  birefringent  opticd 
wedge  and  polarizing  filter.  He  foimd  that,  using  a  CW  laser  source  and  a  12  bit 
resolution  image  sensor  256  pixels  long,  the  minimum  detectable  phase  shift  was  Xy4000 
and  the  accuracy  was  ~X/2000. 

As  mentioned  previously,  the  image  sensor  used  in  the  experiments  presented  here  is 
a  16  bit  resolution  1316x1012  CCD.  With  10  complete  fringes  imaged  onto  the  sensor, 
there  would  be  131.6  pixels  per  fringe.  This  base  resolution,  in  terms  of  pixels/fiinge,  is 
a  factor  of  five  better  than  in  the  Nakadate  experiments.  No  averaging  techniques  were 
used  by  Nakadate. 

Gathering  the  data  for  this  experiment,  50  images  were  saved  with  the  translation 
stage  at  its  starting  position.  The  translation  stage  was  then  stepped  200nm  and  another 
50  images  would  be  acquired.  The  step  and  acquire  was  repeated  several  times. 

After  the  images  were  captured  for  each  position  of  the  translation  stage,  they  were 
processed  using  die  computer  to  analyze  the  data.  The  analysis  program  reports  the 
phase  change  of  each  image  in  relation  to  a  reference  image,  which  in  each  case  is  simply 
the  first  image  acquired.  Plotted  in  Figure  9  is  each  phase  shift  data  point.  Each  point 
represents  the  phase  of  the  image  relative  to  the  first  image  acquired,  sample  number  0. 
The  plot  is  broken  up  into  sets  of  50  images.  Listed  in  Table  2  are  the  mean  values  of 
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each  of  these  sets  of  50  data  points,  their  standard  deviations,  and  the  differences  between 
the  means. 


Figure  9:Phase  shift  data  fi*om  the  verification  experiment.  Each  set  of  50  data  points 
were  taken  with  the  high-resolution  translation  stage  at  a  fixed  position.  The 
stage  was  stepped  200nm  between  each  set  of  data 


Table  (2).  Values  of  the  mean,  deviation,  and  mean  difference  for  each  data 


set  presented  in  Figure  (9) 


Data  Set 

Mean 

Standard  Deviation 

Difference  of  Meap^.;  i 

1 

0.019 

0.0110 

2 

0.055 

0.0137 

0.036 

3 

0.084 

0.0161 

0.029 

4 

0.116 

0.032 

5 

0.144 

0.0130 

0.028 

6 

0.175 

0.0153 

0.031 

The  mean  values  from  Table  (2)  are  plotted  in  the  following  figure.  The  slope  of  the 
line  fit  (and  the  average  of  the  mean  differences)  reveals  a  step  size  of  0.031  radians 
which  is  roughly  X/200.  This  compares  well  with  the  phase  change  estimated  from  the 
measured  fringe  spacing  of  0.036  radians,  a  difference  of  only  16%.  The  linearity  of  the 
measurement  is  readily  apparent  from  Figure  10. 


16 


Figure  10:  Plot  of  mean  values  of  each  data  set  of  50  points  from  Figure  (9).  A  linear 
fit  to  the  data  is  also  shown 


In  order  to  determine  the  minimum  detectable  phase  change,  a  common  statistical 
analysis  was  performed  on  the  data,  the  T-test.  The  primary  purpose  of  the  T-test  is  to 
determine  whether  the  means  of  two  groups  of  data  differ  to  a  statistically  significant 
degree.  Two  data  sets  with  different  means  may  come  from  different  groups  of  data,  or 
the  observed  differences  may  have  occurred  by  chance  or  sampling  error.*®  The  T-test 
helps  determine  which  is  the  case. 

A  T-test  is  a  test  is,  based  on  the  means,  standard  deviations,  and  sample  sizes  of  two 
data  sets,  returns  a  value  which  indicates  the  probability  that  the  two  samples  are  likely  to 
have  come  from  the  same  underlying  population^®.  The  function  used  in  this  evaluation 
returns  a  value  of  1  if  the  two  data  sets  are  identical,  and  returns  a  value  of  0  if  the  data 
set  means  are  widely  separated  with  small  standard  deviations.  Values  between  0  and  1 
are  returned  for  data  sets  whose  means  are  close  or  that  have  large  deviations. 

A  two-tailed  T-test  analysis  assuming  homogeneous  variances  was  performed  on  two 
sets  of  50  data  points  from  Figure  9.  The  analysis  was  done  in  a  spreadsheet  application 
on  a  personal  computer.  The  minimum  detectable  phase  change  was  estimated  by  finding 
the  separation  of  the  means  at  which  the  T-test  returned  a  value  of  0.05.  This  indicates 
that  there  is  a  95%  probability  that  the  two  samples  come  from  different  populations,  or  a 
5%  probability  that  they  come  from  the  same  population.  The  mean  of  one  sample  was 


■V 


ps 


.r^ 

F7 


R 


b'i 


brought  closer  to  the  mean  of  the  other  sample  through  addition  of  a  constant  to  each  data 
point  in  the  set. 

These  tests  reveal  that  the  95%  confidence  level  is  reached  at  a  mean  separation  of 
.0062  radians,  or  X/1021.  So,  with  a  sample  size  of  50  data  points  for  each  set,  the 
resolution  of  the  system  is  about  X/1000. 

4.  Experimental  Results 

i'rom  the  phase  shift  interferometry  measurement,  groups  of-50  data  points  were 
record'*d  for  the  biased  and  unbiased  conditions.  The  input  biases  range  from  2  to  12V.  In 
the  plot  below  (Figure  11),  data  points  are  presented  for  the  unbiased  (OV)  state,  die  solid 
circles,  and  the  biased  (3V)  state,  the  open  circles.  The  input  intensity  is  1.4x10 ’  W/cm^, 
nearly  two  orders  of  magnitude  higher  than  the  measured  traditional  saturation  intensity 
of  3xl0^wW  (for  a  2.5V  bias). 

A  T-test  analysis  was  performed  on  adjacent  data  groups.  That  is,  groups  one  and 
two  were  tested,  two  and  three,  etc..  The  results  of  this  analysis  are  presented  in  Figure 
12.  Except  for  the  tests  14  and  17  (for  data  groups  14/15  and  17/18),  which  do  not 
indicate  a  significant  difference,  other  seven  tests  show  an  80%  or  higher  confidence 
level.  The  average  of  all  the  (1 -T-test)  values  is  0.62.  " 
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Figure  11:  Phase  shift  data  for  OV  (closed  circle)  and  3V  (open  circle)  bias.  Data  was 
collected  in-groups  of  50  for  each  bias  condition,  with  13  sets  of  OV  and  12  sets 
of  3  V  bias  data  being  collected 
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Figure  12:  Mean  for  the  groups  of  50  data  points  presented  in  the  previous  figure.  The 
solid  bars  represent  the  mean  values  for  the  biased  states  (3V)  while  the 
striped  bars  represent  the  \mbiased  state  (OV)  ,  •  - 


As  a  further  analysis  of  the  data,  the  OV  means  were  separated  from  the  3V  means'. 
The  mean  of  the  means  was  then  found  for  both  bias  states.  This  result  is  most  simply 
presented  (Figure  14)  as  two  data  points,  with  error  bars  representing  the  standard 
deviation  of  the  mean  of  the  means  calculation.  It  is  evident  that  there  is  indeed  a 
measured  shift  in  the  phase  of  the  interference  fringe  pattern.  The  difference  in  these 
values  is  fmmd  to  be  -0.0025  rad. 
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Figure  13:  Results  of  a  T-Test  analysis  on  adjacent  groups  of  phase  change  data.  The 
value  of  (1-T-test)  represents  the  probability  that  the  two  sets  of  data  points 
come  from  different  populations 


19 


1::^ 


i 


fV: 

IS: 


rv-T' 


P 


C'' 
!  ■  ; 


3V 


OV 


Figure  14:  The  mean  of  the  means  was  computed  for  both  the  unbiased  and  biased  state. 

Error  bars  represent  the  standard  deviation.  Gaussian  curves  representing 
normal  distributions  are  included  as  an  aid  to  the  eye 

The  change  in  the  refractive  index  can  be  calculated  from  die  change  in  phase 
using  by 

A«i  =  2;r-^  .  (8) 

where  L  is  the  optical  path  length  in  the  material  thickness  and  X.  is  the  wavelength  of 
light.  For  a  reflective  modulator,  the  optical  path  length  is  double  of  the  material 
thickness.  From  the  measured  phase  shift  of -0.0025  rad,  the  change  in  index  is  then  - 
4.7x10'^,  using  Z,=0.9pm  and  X=1.064pm, 

Compared  to  the  calculated  change  in  refiactive  index  using  Kramers-Kronig  relation. 
Figure  15  (a)  and  (b),  the  measured  change  in  refractive  index  is  an  order  of  magnitude 
smaller.  One  possible  cause  is  the  electric  field  screening  in  the  modulator  under  high 
illumination.  The  data  collected  in  the  fringe-shift  measurement  for  voltages  hi^er  th^ 
3V  do  not  show  a  statistically  significant  signal,  within  the  experimeirtal  resolution.  This 


(a)  (b) 


Figure  15:  (a)  Absorption  and  (b)  index  change  for  InGaAs/GaAs  MQW  modulators 
under  low  intensity  illumination  at  I064nm 
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is  an  indication  of  electric  field  screening.  Unfortunately,  due  to  the  lack  of  tunable  pulse 
laser  of  high  power,  the  detailed  dynamic  behavior  of  the  absorption  and  index  of 
refraction  could  not  be  investigated. 

5.  DX  Center  materials  •  •  mxT 

In  this  section,  we  report  on  our  attempt  to  obtain  refiiactive  index  modulation  m  lU-V 

compound  semiconductor  thin  films  that  incorporated  DX  centers. 

DX  centers  are  point  defects  observed  in  n-type  doped  IH-V  compound  semiconductors. 
They  significantly  reduce  the  carrier  concentration  and  cause  persistent 
photoconductivity  (PPC)  at  low  temperature  [21,22].  Figure  16  shows  the  two  stable 
configurations  of  doping  atom,  e.g.  Si,  in  m-V  compound  semiconductors.  In  A^Gai. 
xAs  DX  state  is  the  ground  state  of  donor  atom  when  x  ^0.22.  Depending  on  the 
composition  and  the  sample  temperature  the  donor  atom  in  its  normal  substitutional  site 
(Figurel6-a)  relaxes  into  the  DX  center  configuration  (Figurel6-b)  by  capturing  two 
electrons  from  conduction  band  (CB).  As  a  result  of  this  relaxation  a  bond  between 
doping  atom  and  one  of  its  As  neighbors  is  broken. 


Figure  16:  Two  stable  configuration  of  Si  atom  in  IH-V  compouiids.  (a)  Si  atom  in 
substitutional  site,  and  (b)  Si  atom  in  DX  center  conflation.  Si  atom 
captures  two  electrons  from  conduction  band  when  it  moves  to  DX  center 
configuration 

The  corresponding  energy  of  these  two  states  can  be  represented  in  a 
configuration  diagram  as  shown  in  Figure  17.  In  a  configuration  diagram  total  energy, 
electronics  and  lattice,  of  doping  atom  versus  configuration  coordinate  is  plotted.  In  the 
diagram  shown  in  Figure  17  the  left  curve  represent  the  nonnal  substitutional 
configuration  and  right  curve  represents  the  DX  center  state,  respectively.  The  thermal 
capture  energy.  Heap,  is  important  to  understand  the  PPC  effect  at  low  temperature.  It 
prevents  equilibrium  being  established  below  a  critical  temperature.  When  the  electrons 
are  optically  excited  from  DX  centers  into  the  CB  at  low  temperature  they  stay  in  CB 
because  the  electrons  do  not  have  enough  thermal  energy  to  overcome  the  capture  bamer, 
hence  PPC  effect.  The  persistence  time  of  the  electrons  in  CB  depends  exponentially  on 
the  temperature  and  the  height  of  the  capture  barrier. 


Figure  17:  Configuration  diagram  showing  total  energy  (electronics  and  lattice)  versus 
configuration  coordinate.  When  Si  atom  captures  two  electrons  fi-om 
conduction  band  (CB)  the  lattice  relaxation  takes  place  and  Si  atom  moves  to  a 
distorted  position  Qdx- 

As  a  consequence  of  the  PPC  effect  the  carrier  concentration  in  CB  increases  in 
the  order  of  magnitude.  This  large  change  in  carrier  concentration  affects  the  optical 
properties  of  the  material.  The  expected  index  change  due  to  the  increase  in  fi-ee  carrier 
concentration  can  be  estimated  by  means  of  plasma  effect.  The  dielectric  constant,  e(®), 
of  a  semiconductor  is  given  by 


£(o))  =  e^ico)  - 


(9) 


where  co  is  the  read-out  beam  fi-equency,  and  (Op  is  the  plasma  resonance  firequency  and 


given  by 
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(10) 


Here  m*  is  the  effective  mass  of  electron  in  CB,  and  N  is  the  fi'ee  carrier 
concentration.  The  expression  for  the  change  in  refractive  index  as  a  function  of  the  net 
change  in  carrier  concentration  AN  can  be  expressed  as  follows  [3]. 

(11) 

n^m  CO 


where  n©  is  the  refractive  index  in  absence  of  fi'ee  carriers. 

In  this  work  we  studied  DX  centers  in  a  Si  doped  Alo.3oGao.7oAs  material  grown  by 
molecular  beam  epitaxy  (MBE)  for  a  potential  use  in  index  modulation.  In  this  material 
large  persistent  photoconductivity  (AN=2.1xlO’®  cm*^)  is  observed.  We  measured  the 
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anisotropy  of  the  sample  resistivity  (Ri/R//  >  10'*)  resulted  from  modulation  of  the  free 
carrier  density  by  optical  excitation  in  striped  pattern. 

The  sample  studied  in  this  work  is  a  1  pm  thick  Si  doped  AlojoGao.roAs  grown  by 
molecular  beam  epitaxy  (MBE)  on  a  semi-insulating  GaAs  substrate.  An  additional  1  pm 
undoped  Alo.30Gao.70As  layer  is  also  grown  between  substrate  and  doped  Al0.30Ga0.70As 
layer  to  prevent  quantum  well  formation  between  GaAlAs  and  GaAs  interface.  After  the 
growth  the  sample  is  cut  into  a  5x5  mm  sample  size  to  do  electrical  measurement  in  a  van 
der  P'iaw  geometry.  Ohmic  contacts  are  made  to  the  sample  with  In.  - 
Figure  18  shows  the  measured  PPG  data.  In  this  measurement  the  sample  is  cooled  down 
from  room  temperature  to  25  °K  in  the  dark.  At  25  “K  it  is  uniformly  illuminated  with  a 
white  light  source  20  minutes  to  excite  all  the  DX  centers.  After  the  light  is  turned  off, 
we  waited  an  additional  10  minutes  to  make  sure  that  the  sample  has  reached  a  steady 
state.  The  carrier  concentration  is  then  remeasured.  The  net  change  in  earner 
concentration  below  77  ‘X  is  around  AN=2.1xlO’*  cm*^  This  increment  in  c^er 
concentration  is  persistent  as  long  as  the  temperature  is  kept  below  the  critical 
temperature,  which  is  120  °K  for  GaAlAs  based  materials. 


10^® 
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Figure  18:  The  persistent  photoconductivity  (PPG)  curve  for  Alo.30Gao.70As: Si.  The 
sample  is  cooled  down  from  room  temperature  to  25  ‘TC.  At  25  K  it  is 
uniformly  illuminated  20  minutes  with  a  white  light  source  to  excite  all  the  DX 
centers. 

As  mentioned  above  the  modulation  of  the  carrier  density  leads  to  large  refractive  index 
change  as  a  result  of  the  plasma  effect.  The  local  nature  of  the  DX  centers  can  be 
exploited  by  illuminating  the  sample  with  a  striped  pattern  parallel  to  one  of  the  sample 
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edges.  The  illuminated  portion  of  the  sample  will  have  larger  carrier  concentration  than 
that  of  dark  region.  The  difference  in  the  refiactive  index  due  to  this  spatial  charge 
modulation  between  these  two  regions  can  be  estimated  by  using  eqs.  (1 1). 


Figure  19  shows  the  experimental  setup  used  to  measure  such  index  change.  The  sample 
is  placed  in  an  optical  cryostat  which  allows  us  to  do  low  temperature  optical  and 
electrical  measurements.  In  this  experiment  a  striped  pattern  in  different  grating  periods 
is  placed  in  front  of  the  sample.  The  sample  is  illuminated  with  a  HeNe  laser  (X=6328  A) 
which  has  enough  optical  energy  to  excite  electrons  from  DX  centers.  This  ojptical 
illumination  in  a  striped  pattern  creates  modulation  of  the  free  carriers  along  the  side  of 
the  sample.  The  sample  conductivity  is  measured  along  perpendicular  (Ri)  and  parallel 
(R// )  directions  using  a  van  der  Pauw  geometry. 
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Figure  19:  Experimental  setup  used  to  measure  the  modulation  of  the  carrier  density.  The 
sample  is  illuminated  with  a  HeNe  laser  using  a  striped  pattern  parallel  to  one 
of  the  sample  edges. 


The  measured  resistivity  in  perpendicular  and  parallel  directions  as  ^function  of  grating 
space  is  shown  in  Figure  20.  This  figure  also  shows  the  magnitude  of  the  resistivity 
anisotropy.  When  the  sample  is  uniformly  illuminated  no  anisotropy  in  conductivity  is 
observed.  However,  if  the  sample  is  illuminated  with  a  striped  pattern  parallel  to  one  of 
the  sample  edges  a  large  anisotropy  in  resistivity,  up  to  lO'*,  is  observed.  It  is  possible  to 
measure  the  effect  of  this  charge  modulation  optically  in  our  setup  with  a  slight 
modification.  Linke  et  al.  optically  measured  a  refractive  index  change  An=1.2xl0‘  in  an 
Alo.29Gao.3iAs:Si  sample  doped  to  a  value  of  4.0x10**  cm*^  [24]. 
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Figure  20:  Modulation  of  the  resistivity  in  Alo.30Gao.70As  as  a  function  of  grating  period. 

There  is  no  anisotropy  in  resistivity  when  the  sample  is  uniformly  illuminated. 
The  anisotropy  in  resistivity  up  to  K±fRji  >10^  is  measured  for  large  grating 
size. 

6.  Conclusions 

We  have  demonstrated  experimentally  the  refractive  index  modulation  in  MQW 
materials  beyond  the  absorption  saturation.  The  measur^  index  change  at  high  optical 
illumination  is  an  order  of  magnitude  smaller  than  estimated  from  the  low-intensity 
absorption  coefficient,  which  may  indicate  the  presence  of  electric  field  screening  under 
the  high-intensity  illumination. 

As  discussed  in  Section  2,  the  electric  field  screening  is  caused  by  the  imbalanced 
escaped  rates  of  electrons  and  holes  in  the  wells.  In  general  electrons  are  lighter  than 
holes  in  any  material,  the  screening  is  more  likely  caused  by  holes  left  in  the  wells. 
Therefore,  the  materials  with  higher  potential  barrier  for  electrons,  such  as  GaAs  related 
component  materials,  are  better  suited  for  the  devices  operating  at  high  intensity.  It  has 
been  reported  that  the  incorporation  of  Nitrides  in  ni-V  semiconductor  materials  results 
in  the  increasing  of  conduction-band  offsets  in  the  MQW  materials.  Further  research 
could  be  conducted  in  Nitride  related  materials  if  the  MQW  devices  are  to  be  used  to 
modulate  phase  under  high-intensity  illumination. 

Finally,  the  phase-shift  interferometry  method  developed  during  the  course  of  this 
work  is  proved  to  be  a  high-resolution  technique  to  determine  the  index  change  in  a 
material.  It  has  the  resolution  of  ^72500. 

Due  to  the  local  nature  and  PPG  effect  the  DX  centers  based  index  modulation  could 
be  important  for  optical  devices.  The  refraction  index  of  the  material  that  contains  DX 
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centers  can  in  principle  be  modulated  by  excitation  of  these  defects,  hi  this  work  we 
measured  the  spatial  modulation  of  the  carrier  concentration  in  Al0.30Ga0.70As  material 
and  observed  a  large  anisotropy  in  conductivity  (Rx/R//>  10^)  at  77K. 

Although  we  have  tried  to  measure  the  refractive  index  of  our  materials  with  the 
above-describ^  set-up,  we  were  not  ^le  to  demonstrate  any  optical  modulation  in  our 
DX  materials.  We  suspect  that  due  to  the  small  thickness  of  our  films  the  phase 
retardation  was  too  small  for  our  set-up  to  sense.  Further  studies  are  needed  in  this  area. 
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